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Abstract 

A fully numerical method to calculate loop integrals, a numerical 
contour-integration method, is proposed. Loop integrals can be inter¬ 
preted as a contour integral in a complex plane for an integrand with 
multi-poles in the plane. Stable and efficient numerical integrations 
an along appropriate contour can be performed for tensor integrals as 
well as for scalar ones appearing in loop calculations of the standard 
model. Examples of 3- and 4-point diagrams in 1-loop integrals and 
2- and 3-point diagrams in 2-loop integrals with arbitrary masses are 
shown. 

Moreover it is shown that numerical evaluations of the Hyperge¬ 
ometric function, which often appears in the loop integrals, can be 
performed using the numerical contour-integration method. 
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1 Introduction 


In future collider experiments such as at the LHC and ILC, the standard 
model will be checked with very high precision and a signal of new physics 
will be searched for through a tiny difference between experimental mea¬ 
surements and theoretical predictions. The theoretical uncertainty must be, 
at least, one order of magnitude smaller than the experimental one. The 
theoretical prediction is obtained based on the perturbative calculation of 
the quantum field theories. The theoretical calculation of a higher order 
correction with one- or two-loops must be performed in order to meet with 
experimental requirements. A loop integration is one of the critical issues 
of a computation of these higher order corrections. Basically there are two 
methods to perform the loop integration, an analytical method and a nu¬ 
merical one. Though the analytical method gives, in principle, fast and 
stable results, for many cases it is hard to give a compact expression for the 
one-loop integral and it is very difficult to perform two-loop computations 
for the standard model with multiple energy scales. 

Numerical methods have also been investigated to perform loop inte¬ 
grals. The symmetrical sampling method was proposed in ^ by Oyanagi et 
al. in 1988 and has been investigated through a series of papers [3 El 13 • An¬ 
other numerical method, the hybrid method jK], has been proposed. It is to 
perform a part of the multi-dimensional integrations analytically up to the 
integrand being a logarithmic form. A remaining integration was performed 
numerically by Monte Carlo method. Both methods are tuned to calculate 
up to two-loop/three-point functions with arbitrary masses [H]- 

Recently another numerical methods using the Sato-Bernstein-Tkachov 
(SBT) relation ^ is proposed and utilized intensivelv |l()[ ITTl IT^ - A recent 
review of the method using the SBT relation can be found in m- 

Yet another numerical method has been proposed by de Doncker|14j. 
an ‘e-algorithm’. This method takes the infinitesimal imaginary parameter 
appearing in a denominator of the loop integral at finite values, and takes its 
extrapolation to zero. It is demonstrated that this method can give precise 
results for a non-scalar one-loop/three-point digram^3- 

We propose a new method of loop integrations such as a ‘numerical 
contour-integration (NCI) method’ in this report. The loop integral can 
be interpreted as a contour integral in a complex plane for integrands with 
multi-poles in the plane. We show that the stable and efficient numerical 
integrations along an appropriate contour for tensor integrals as well as 
scalar ones appear in loop calculations of the standard model. This method 
is applicable, so far, up to two-loop/three-point functions with arbitrary 
masses and arbitrary polynomials of Feynman parameters in the numerator 

^This relation was introduced at first by Sato[2| relating to prehomogeneous vector 
spaces and was generally proved by Bernsteinj^]. An application of this relation to the 
loop integral was proposed by Tkachov|3. 
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of the integrand. The basic idea of this method and some numerical results 
are shown in this report. 


2 Feynman parameter representation of tensor in¬ 
tegrals 


The tensor integral of a massless one-loop N-point graph with rank M < N 
in a space-time dimension of n = 4 — 2e can be written as, 


rp{N) 

^ fi - ■ ■ ly 


M 


(Pk k,. 


(27r)"z A 1 A 2 ■ ■ ■ An ’ 


where 


Ai — {k — Tnf iO, 

i 

Si = ^Pj, so = 0, 


and Pi is a four momentum of an i’th external particle (incoming), k^ a loop 
momentum, and m* the internal masses. An infinitesimal imaginary part 
(iO) is included to obtain analyticity of the integral The momentum 

integration can be done using Feynman’s parameterization which combines 
the propagators. After the momentum integration, an ultra-violet pole is 
subtracted under some renormalization scheme. Finally the tensor integral 
can be expressed using integrals of the type 


where 



x\ 




{Dn - iO) 


N-2 ’ 
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N-l 
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Xi 


dXN-l- 
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An explicit form of the numerator, and a relation between and can 
be obtained straightforwardly depending on the diagram^^l- The remaining 
task is to perform the parametric integration of in Eq. o, which is called 
a tensor integral in this report. (When the numerator of the integrand is 
unity, it is called a sealar integral.) 


3 one-loop/three-point 


A general formula of the tensor integral for the one-loop/three-point function 
with arbitrary masses is^ 


/( 3 ) _ 


/o Jo Do — *0 
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Figure 1: l-loop/3-point diagram 


D 3 = (Mix — M 2 y)‘^ — rxy + m\ix + m‘^ 2 y + ^ 2 ; 

where 


r = s-(Mi+M 2 ) 2 , 

2 7\ /r‘2 I 2 2 

= —Ml + nil—m2, 

2 7 \ //'2 I 2 2 

77132 = -M2 +7773 - 7772, 

Mf = 

M| = pI 

S = ipi+P2f- 


A numerator of the integrand, f{x,y), can be any polynomial of Feynman 
parameters x and y with rank M < 3. Momentum and mass assignments 
are shown in Figure 1. 

The integration region is a 2-dimensional simplex with three sides. Those 
sides are 


Lx = {(x,7/)|0 < X < 1,7/ = 0}, 

Ly = {(x,7/)|x = 0,0 < y < 1}, 

Lxy = {{x,y)\y = l - x,Q <x <l}. 

Singular points (second pole) of the integrand form a parabola (or hyper¬ 
bola), which is 


'P = {(x,y)|T)3 = 0}. 
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The focal point of the parabola is on the opposite side of the origin with 
respect to the parabola line. When one integration-variable is fixed, the 
remaining one-parameter integration may intersect a pole on the line V. 
This integration can be done numerically as a contour integral along an 
appropriate contour with correct analyticity. For a stable integration, it 
must be avoided that the contour intersects more than one pole, or that a 
pole is at the end point of the contour. These requirements are satisfied by 
taking an appropriate coordinate system instead of a simple (x, y) coordinate 
of Feynman parameters. Basically we take polar-coordinates with an origin 
on the side L^y Coordinate systems we used are categorized by the following 
cases; 

1 . M[{L^iJLyiJL^y)r\V] = 0 : 

AA[^] are number of elements of a set A. There is no singular points 
in the integration region. The origin is set at the nearest point to the 
line V on the side L^y 

2. M[L^y r\V]=l ^ M[{L^ U Ly) n T’] = 0: 

The origin of the coordinate system is set at a point such as {Lxyf^V}. 
Though a pole is at the origin, the singularity is canceled out against 
the Jacobian on the numerator, i.e. the measure at the origin is zero in 
the polar-coordinates. Then the integration is free from singularities. 

3. M[Lxy r\V]=2 ^ M[{Lx U Ly) n T’] = 0: 

The origin is set at the center of two elements of {Lxy H "P}. Then 
the contour along r intersects the line V only once for any value of 
(j). We can take an appropriate contour avoiding the pole with correct 
analyticity. 

4. J\f[Lxy n P] = 2 A N[Lx n P] = 2 A N[Ly n P] = 2: 

The integration region is divided into three regions by two lines, one 
connects between centers of two elements of Lxy H P and those of 
Lx n P, the other connects between centers of two elements of Lxy H P 
and LyDV. The origins of the coordinate frames are set at the three 
corner of the simplex, (0,0), (1,0) and (0,1). In each region, the 
contour along r intersects the line P only once for any value of (j). 
We can take an appropriate contour avoiding the pole with correct 
analyticity. 

Numerical integrations are performed based on the ‘Good Lattice Point 
(GLP) Methoddni’, which uses a deterministic series of numbers. For 
smooth functions it shows very efficient convergence compared with a Monte 
Carlo integration based on random sampling. Moreover both imaginary and 
real parts can be integrated simultaneously. 

Here we show several examples of one-loop/three-point functions. 
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Case l:mi = m 3 = Mi = M 2 = 150 GeV, m 2 = 91 GeV, f{x,y) = 1. 


Table 1 



real/imag. 

analytic [17j 

NCI result 

error 

calls 

310 

real 

imag. 

0.1110 X 10“^ 
0.7162 X 10“^ 

0.110975 X 10“^ 
0.716174 X 10““^ 

0.23 X 10-^'J 

350 

500 

real 

imag. 

0.2285 X 10“^ 
0.4731 X 10“^ 

0.228498 X lO”'^ 
0.473072 X 10“^ 

0.39 X 10“^'^ 

350 

1000 

real 

imag. 

-0.6103 X 10“^ 
0.1551 X 10“^ 

-0.610293 X 10“^ 
0.155066 X lO-'^ 

0.22 X 10“^^ 

582 


This is an infrared-divergence free integral with only heavy particles 
involved. Very good agreements can be obtained between analytical 
results (formulae given in ref.|17|i and the NCI method as shown in 
the above table with only several hundred sampling points. 

• Case 2:mi = m 3 = Mi = M 2 = mg, ^/s = 1000 GeV, f{x,y) = 1. 


Table 2 


m2{GeV) 

real/imag. 

analytic jT7j 

NCI results 

error 

calls 

10 “^ 

real 

imag. 

-0.641175 X 10'^ 
0.115741 X 10-3 

-0.641175 X 10-3 
0.115740 X 10-3 

0.20 X 10 -y 

6732 

10 -^ 

real 

imag. 

-0.907973 X 10-3 
0.144676 X 10-3 

-0.907973 X 10-3 
0.144674 X 10-3 

0.22 X 10-'^ 

6732 


real 

imag. 

-0.117483 X 10-^ 
0.173611 X 10-3 

-0.117479 X 10-^ 
0.173618 X 10-3 

0.10 X 10-^^ 

28621 


This is an infrared-divergent case with a fictitious photon mass of 
10“^ to 10“® GeV. Very good agreements can be obtained between 
analytical results and the NCI, even for the infrared-divergent case 
with several thousand sampling points. 

• Case 3:mi = m 3 = 80 GeV, Mi = M 2 = 170 GeV, M 2 = 5 GeV, 
Vs = 357 GeV 


Table 3 


f{x,y) 

real/imag. 

FF[IB 1 

NCI results 

error 

calls 

1 

real 

imag. 

-0.623169 X 10-4 
0.132340 X 10-4 

-0.623169 X 10-4 
0.132340 X 10-4 

0.18 X 10-44 

1567 

X 

real 

imag. 

-0.171520 X 10-4 
0.579780 X 10-3 

-0.171520 X 10-4 
0.579780 X 10-3 

0.43 X 10-^ 

1567 

X‘^ 

real 

imag. 

-0.859550 X 10-3 
0.487062 X 10-3 

-0.859550 X 10-3 
0.487063 X 10-3 

0.13 X 10-y 

1567 

X X y 

real 

imag. 

-0.324021 X 10-3 
0.296541 X 10-3 

-0.324021 X 10-3 
0.296540 X 10-3 

0.11 X 10-y 

1567 
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Figure 2: l-loop/4-point diagram 


The tensor integration can be done as well as the scalar one as shown in 
Table 3. NCI results show a very good agreement with those from the FF 
package [T^ developed by van Oldenborgh. 


4 one-loop/four-point 


A general formula for the one-loop/four-point function with arbitrary masses 

is|Hl 


/( 4 ) 

Da 


rl-z 


dz 


dy 




L >4 — iO ’ 


{Mix - M 2 ?/)^ - rxy + mhx -|- 77142 ?/ + 77143 , 


where 


t = {P2+P3f, 

u = {pi+p^f, 

mil = —Mf + 7714 — 7712 — 2;(s + 71 — 2Mf — M| — M|), 

77142 = —M| + 77l| — + M| — M|), 

"1I3 = 7712 -|- (—-Tf| + 7714 — 7712)2; + tz ^ 


A numerator of the integrand, f{x, y, z), can be any polynomial of Feynman 
parameters, x, y and 2 ; with rank M < 4. Momentum and mass assign¬ 
ments are shown in Figure 2. The integration region is a 3-dimensional 
simplex. In order to perform the contour integral very efficiently in this 
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three-dimensional integration-parameter space, we have introduced a new 
coordinate system as shown in Figure 3, named ‘Wedge coordinate system’. 
A relation between an usual Cartesian coordinate system and the wedge 



Figure 3: Parametrization of the wedge coordinate 


coordinate system is given as 


X = r cos 4>' sin { 9 + 


where 


y = r cos (p cos \^9 -|- 
^ = {Xmax - r) sin (j)', 


tan (h' = cos 9 tan ( 


TT 

tt' 


I tan^ 6 1 

+ 


2 ' 2cos2 6»' 

The Feynman parameter integration now becomes 

/(4) = 


/o Jo Jo 


Di — iO 

t-Tmax r7r/4 COS 9 COS^ 4>' f (x, 9, (j)) 

/ dx d9 dcf) x{Xmax - X)- 

10 J-Tr /4 Jo 


COS^ (p D 4 — iO 


When cj) is fixed, the remaining integration region is a triangle. The inter¬ 
section between this triangle and a hyper-surface which satisfies £>4 = 0 
is a parabola (or hyperbola) similar to the one-loop/three-point case. The 
variable r is extended into the complex plane. When both cj) and 9 are fixed, 
we can take an appropriate contour for r which crosses a pole only once for 
any </> and 6 . 

Here we show several examples of one-loop/four-point functions. 

• Case l:e+e- ^ W+W. 

= 80 GeV, m 2 = 0, = 91 GeV, Mi = M 2 = me, 

Ms = M 4 = 80 GeV, ^ = 500 GeV, 9 = L{Pi,Pa). 


Table 4 


cos 9 

real/imag. 

FFdHl 

NCI results 

error 

calls 

-0.5 

real 

imag. 

-0.459774 X 10“^ 
0.711583 X 10-9 

-0.459774 X 10-9 
0.711582 X 10-9 

0.67 X 10-^^ 

1567 

0.0 

real 

imag. 

-0.588538 X lO-^ 
0.962985 X 10-9 

-0.588539 X 10-9 
0.962985 X 10-9 

0.10 X 10-^^ 

1567 

0.5 

real 

imag. 

-0.849170 X 10-9 
0.155695 X 10-® 

-0.849167 X 10-9 
0.155695 X 10-® 

0.10 X 10-^^ 

2553 


Though a massless particle (neutrino) appears in the loop in this case, 
this is infrared-divergence free. The NCI with several thousand sam¬ 
pling points gives very good agreement with FF. 

• Case 2:e'''e“ ^ zz. 

mi = m 3 = m 4 = me, m 2 = 91 GeV, Mi = M 2 = me, M 3 = M 4 = 91 
GeV, Vs = 500 GeV, 9 = 1{Pi,Pa). 


Table 5 


cos 9 

real/imag. 

FT [IB] 

NCI results 

error 

calls 

-0.5 

real 

imag. 

-0.89532 X 10-9 
-0.23060 X 10-9 

-0.89540 X 10-9 
-0.23049 X 10-9 

0.54 X 10-^^ 

2M 

0.0 

real 

imag. 

-0.11548 X 10-« 
-0.33447 X 10-9 

-0.11549 X 10-« 
-0.33433 X 10-9 

0.78 X 10-1^ 

2M 

0.5 

real 

imag. 

-0.17129 X 10-« 
0.60873 X 10-9 

-0.17134 X 10-^ 
0.60862 X 10-9 

0.96 X 10-^^ 

2M 


A light particle (an electron) exists in the loop in this case. Although 
it requires high statistics of about two million sampling points even by 
the GLP, the result shows good agreement with FF. 

Though only scalar integrations are shown here, tensor integrations with an 
arbitrary polynomial in the numerator can be done as well. 
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Figure 4: 2-loop/2-point diagram 


5 two-loop/two-point 


A general formula of the two-loop/two-point function with arbitrary masses 

islH 


/( 2 ) 


where 


dxidx2dx^dxiidx^6 1 



fjxi) 

DC-iO' 


D = -p^ {X5{xi + Xs) {X2 + Xi) + {xi + X2) X3X4 + {X3 + X4,) X1X2) + CM'^, 

C = {xi + X2 + X3 + X4) X5 + (xi + X2) (X3 + X4) , 

5 

Xjmf, 

i=l 


and rrii are the internal masses. A particle number assignment is shown in 
Figure 4. Here f{xi) is any polynomial of Feynman parameters Xi with rank 
M <2. 

Feynman parameter integration with four independent variables must be 
performed after X 5 integration eliminates a delta-function. It is very hard to 
make any intuitive optimization for the integration because four-dimensional 
space is beyond our imagination. We simply employ the following coordinate 
system: At first one-point on the edge-simplex oi xi + X 2 + X'i + x^ = 1 
is chosen. Then take a variable r as a distance between a origin of the 
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coordinate (xi = X 2 = xs = 2:4 = 0 ) and the point on a straight line 
from the origin to the point in the edge-simplex. This r is extended into 
the complex plane and chosen to follow an appropriate contour to avoid 
singularities. Here we show two examples of two-loop/two-point functions. 

• Case l:mi = m 2 = m 3 = m 4 = mo = 150 GeV and ms = 91.17 GeV. 


r 

!'able 6 

/rriQ 

real/imag. 

Kreimer)20j 

NGI results 

2.0 

real 

2.664 

2.6647(2) 

4.0 

real 

25.85 

25.87(2) 

4.1 

real 

21.506 

21.47(3) 


imag. 

12.16 

12.15(3) 

4.3 

real 

15.246 

15.19(3) 


imag. 

17.013 

17.29(3) 

5.0 

real 

3.668 

3.63(4) 


imag. 

19.503 

19.48(4) 

7.0 

real 

-6.7939 

-6.78(4) 


imag. 

13.918 

13.89(4) 

10.0 

real 

-8.793 

-8.792(1) 


imag. 

8.865 

8.846(1) 


In this case, integrations are done by Monte Carlo method using 
BASES [E]. NCI results show good agreement with the analytical 
results by Kreimerj 2 ()j. within the statistical error of the Monte Carlo 
integration. 

• Case 2:mi = 1 GeV, m 2 = 2 GeV, m 3 = 4 GeV, m 4 = 5 GeV, ms = 3 
GeV. 


Table 7 


p\GeV'^) 

real/imag. 

Bauberger et al.pT) 

NCI results 

0.1 

real 

-0.287238(3) 

-0.287240(6) 

0.5 

real 

-0.294592(3) 

-0.294594(6) 

1.0 

real 

-0.304521(3) 

-0.304523(6) 

5.0 

real 

-0.452520(3) 

-0.452527(9) 

10.0 

real 

imag. 

-0.488153(2) 

-0.353217(2) 

-0.48807(5) 

-0.35307(4) 

50.0 

real 

imag. 

0.173901(2) 

-0.118080(2) 

0.17397(2) 

-0.11805(2) 


This is an example with a very simple mass assignment. The result of the 
NGI method is compared with that from the analytical calculation obtained 
by Bauberger and BbhmPT|. and gives also very good agreement. The same 
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Figure 5: 2-loop vertex diagram (non-planer type) 


example is treated using the SBT relation by Passarino and Uccirati|llj with 
very good agreement too. 

Though only scalar integrations are shown here, tensor integrations with 
an arbitrary polynomial in the numerator can be done as well. 


6 two-loop/three-point 


A general formula of the two-loop/three-point function (non-planar config¬ 
uration) with arbitrary masses can be expressed as 

f{yi,Zj) 


= F / dzidz 2 dz ^5 (1 - az) ziZ 2 Z^ [ dyidy 2 dy 3 
o JO J-i 


{Do - iO) 


2 ’ 


after some transformation of the original Feynman parametrizationjn]. Here 
the denominator is 


where 


A = - 


Do = ^^A-y 

^ T • ^ + c, 

( yi\ 


It = 

CO to 


^ -zl{z2 + Zo)si 

ZlZ2Zo{-Sl - S 2 + So )/2 

ZlZ2Zo{-Si + S 2 - So )/2 

ZlZ 2 Zo{-Sl - S 2 + So )/2 

-zl{zo + Zi)s 2 

ZiZ2Zo{+Si - S 2 - So )/2 

V ZlZ 2 Zo{-Sl + S 2 - So )/2 

- S2 - So )/2 

-zl{zi -h Z 2 )so 
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c 


h 


^ zi{m\ — ml) 

Z2lml - ml) 

V Z 3 {ml - mj) 



U 



Zisi + Z2S2 + 


Z1Z2 + 5:2^3 + 


2 (m| + ml)zi - 2 {ml + ml)z 2 - 2 {m\ + ml)z 3 , 


and Si = Pi, and mj are the internal masses. Here gives a transposed 
vector of if. For the y integration, we employ the SET relation as 


Dl 


1 

P 


1+2 -^)- dy 


1 
J 


P = c-{^lPA-^lPy 


V = 
= 


-a-^T 


/ d d d 
\dyi' dy2' dys 


The yi are left as real variables. In order to avoid singularities in /3“^ and 
the Zi are extended into the complex plane. After the Z 3 integration to 
eliminate the 5-function, the polar coordinates, zi = sin^^, Z 2 = cos 9z 
are introduced. This is extended into the complex plane and chosen along 
an appropriate contour to avoid singularities. Here we show one example of 
a two-loop/three-point function. 


• mi = m 2 = m/i = ms = 150 GeV, m 3 = mg = 91.17 GeV, and 
si = s 2 = 1502 GeV 2 . 


Table 8 


S 3 /mi 

Fujimoto et al.0 

NCI results 

4.5 

2.09(2) X lO^y 

2.049(2) X 10“^ 

5 

1.43(1) X 10“^ 

1.464(2) X 10“^ 

12 

-8.37(3) X lO-^'^ 

-8.24(4) X 10“^^^ 

20 

-5.61(3) X lO-^'^ 

-5.60(1) X lO-^'J 


The Monte Carlo integration package, BASES, is used for the numerical 
integration. The results agree well with those in refs. 0122] 

Though only scalar integrations are shown here, tensor integrations with 
arbitrary polynomial in the numerator can be done as well. 

7 Numerical evaluations of the Hypergeometric 
functions 

So far we discussed in previous sections loop-integrals in the standard model 
with arbitrary masses. In this section we would like to treat a massless 
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theory such as QCD. A general one-loop/four-point function in a massless 
theory can be expressed by the following tensor integrals: 

>{s,t-,no„ny,nz) = (^dvr/r j dx dy dz ’ 

Di = —xzs — y{l — X — y — z)t, 


where s = {pi +P 2 )‘^, and t = {pi +^ 4 )^. Here we set all particles massless. 
In order to cure an infrared divergence, the space-time dimension is set to 
be n = 4 -|- 2sm after the MS renormalization. The tensor integration can 
be done analytically and be represented by a finite number of terms with 
Beta and Hypergeometric functions asjSnj; 


X 

X 

-h 

X 

s = 

i = 

t = 
u = 


t( 4) r , \ 1 B{nx + EiRTTiy + riz + Sjji) 

T >{s,t-,nx,ny,nz) = —-^^- 

I - 


s t 


Attp? 


Sir 


s 


Tlx 


Wj=i{nx- j + eiR) 


B{1 -h riz^Ux + ny + Sir) 


u 


2Fi\l + Ux, nx + nz + Sir, l + nx + ny + nz + sir,-- 


— Z\ ^IR / _ 


Airp^ 


1=0 




i ni 


- j - nx) 
^) n7'=i(^ — 3 +^ ir ) 


B{1 -|- Uy, I + nz + Sir) 


2 F 1 + l,l + riz + Sir, 1 + I + riy + Uz + sir, — 

s + iO, 
t + iO, 
t — iO, 

u + iO = {pi + ^0. 


u 


For the infrared-finite case, evaluation of the loop integral can be done if we 
can perform the numerical calculation of the Hypergeometric function of a 
type 


2 T 1 (1 -I- /, / -I- m, 1 -I- / -I- m -I- n, z) 


r{l +I+ m + n) F 

T{1 + m)T{l + n) Jo ^ {1 - zt)^F 


where l,m,n > 0 are integer numbers and z a complex variable. This 
integration can be performed easily using the NCI. 
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1 

m 

n 

real/imag. 

M athematica |Uj 

NCI 

calls 

1 

1 

1 

real 

imag. 

-0.1453322029 x 10“^ 
-0.1507964474 x 10° 

-0.1453322029 x 10“^ 
-0.1507964474 x 10° 

4149 

1 

2 

3 

real 

imag. 

0.8417767168 x 10-^ 
-0.2290221044 x 10° 

0.8417767169 x 10-^ 
-0.2290221044 x 10° 

4149 

2 

1 

1 

real 

imag. 

0.1087664688 x 10“^ 
0.2638937829 x lO”! 

0.1087664688 x 10“^ 
0.2638937829 x lO”! 

4149 

2 

3 

4 

real 

imag. 

-0.2890568082 x 10“^ 
0.5578978464 x 10-^ 

-0.2890568082 x 10“^ 
0.5578978464 x 10“^ 

6732 

3 

1 

2 

real 

imag. 

-0.1026721798 x 10“^ 
-0.5654866776 x lO-^ 

-0.1026721798 x 10“^ 
-0.5654866776 x lO-^ 

6732 

3 

4 

5 

real 

imag. 

0.8121358824 x 10“^ 
-0.1274636264 x 10“^ 

0.8121358824 x 10“^ 
-0.1274636264 x 10“^ 

6732 


Table 9 numerical calculation of the Hypergeometric function at 2 = 10. 


The results from the NCI method with the CLP numerical integration at 
z = 10 + iO are shown in Table 9 comparing them with those obtained 
by Mathematical^. The ten-digit agreement with Mathematica can be 
obtained with only four to seven thousand sampling points in the CLP. 

8 Summary 

We have developed a fully numerical method, named ‘numerical contour 
integration (NCI)’ method, to calculate loop integrals. Loop integrals can 
be interpreted as a contour integral in a complex plane for an integrand with 
multi-poles in the plane. For the one-loop/three-point case with arbitrary 
masses in the standard model, both scalar and tensor integrals have been 
obtained by the NCI method with five to seven digits accuracy. For the one- 
loop/four point case, scalar integrals have been calculated with about 0.1% 
accuracy. For the above two cases, the ‘Good Lattice Point’ method is used 
for the numerical integration. Those results are compared with analytical 
calculation or FF and show a very good agreement. The NCI method has 
been applied also to two-loop integrals. For two- and three-point integrals 
at the two-loop level, it is demonstrated that the NCI method can give 
numerical results of the loop integrals with good accuracy using the Monte 
Carlo integration package BASES. Those results show a good agreement 
with previous calculations. The extension to the general tensor cases is 
straight forward, since the NCI is a purely numerical method. Moreover 
it is shown that the numerical evaluation of the Hyper geometric function, 
which appears in the one-loop/four-point tensor-integrals in massless QCD, 
has been performed with ten-digit accuracy using the CLP method. 
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